Quantum fluctuations in the image of a Bose gas 
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We analyze the information content of density profiles for an ultracold Bose gas of atoms and 
extract resolution limits for observables contained in these images. Our starting point is density 
correlations that we compute within the Bogoliubov approximation, taking into account quantum 
and thermal fluctuations beyond mean-field theory. This provides an approximate way to construct 
the joint counting statistics of atoms in an array of pixels covering the gas. We derive the Fisher 
information of an image and the associated Cramer-Rao sensitivity bound for measuring observables 
contained in the image. We elaborate on our recent study on position measurements of a dark soliton 
[Negretti et al, Phys. Rev. A 77, 043606 (2008)] where a sensitivity scaling with the atomic density as 
n~'^' ■* was found. We discuss here a wider class of soliton solutions and present a detailed analysis 
of the Bogoliubov excitations and the gapless (Goldstone) excitation modes. These fiuctuations 
around the mean field contribute to the noise in the image, and we show how they can actually 
improve the ability to locate the position of the soliton. 

PACS numbers: 03.75.Nt, 06.20.Dk, 37.25. +k 






en 
> 

p 
o 

OO 

o 



X 



I. INTRODUCTION 

de Broglie waves of massive particles are very sensitive 
to perturbations and may serve as efficient probes for 
electromagnetic fields Ij, earth's rotation [3], Casimir 
forces Q (or in general to detect weak forces Q), and 
particle properties as, for example, the refraction index 
of a buffer gas ^ or the electric polarizability of an atom 
[y|. In interferometry, disturbance of the phase of light 
or matter waves in one arm of the interferometer can be 
measured by a displacement of the interference fringes 
with a sensitivity determined by the fringe wavelength 
and the signal-to-noise ratio (SNR). The SNR is under 
many circumstances given by the standard shot noise, 
leading to a resolution that scales with 1/^/N, where N 
is the number of detected atoms. 

The shot noise limit, however, can be beaten with en- 
tanglement and squeezing [3] , proving that it is not a fun- 
damental limit and in the scenario of Ref.[8i], for exam- 
ple, a particular entangled preparation scheme has been 
shown to give a phase error scaling as l/N^'"^, while in 
principle the Heisenberg limit should provide the opti- 
mum sensitivity with a phase error scaling as 1/N [9|. 
Many analyses of nonclassical metrology with quantum 
objects have dealt with the situation of particles or fields 
that may be prepared in suitably entangled states, and 
entangled measurements may be used after they have ex- 
perienced the interaction of interest [10| . Conversely, en- 
tanglement created by many-body interactions improves 
the precision of estimating the corresponding coupling 
constant, as discussed recently }lll.ll2|. 

In this paper we arrive at a 1/N^^'^ scaling for a dis- 
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placement measurement with particles that are, however, 
not entangled at all. We consider a system of bosonic 
atoms, which are cooled to Bose-Einstein degeneracy and 
prepared in a joint collective quantum state, described by 
mean field theory, i.e., by a Hartree-Fock product state 
wave function. We take for the latter a dark soliton, a 
topological excitation stabilized by atom-atom interac- 
tions. The soliton has a density minimum at a location 
q, and this position can be used to monitor the phase 
acquired in an atomic matter wave interferometer in the 
nonlinear regime [13l . [iJ, |l5[ . A recent experiment has 
demonstrated the relevance of the relative phase in split- 
ting and recombining a Bose-Einstein condensate (BEC), 
although the direct observation of solitons was not pos- 
sible 16]. In a similar way, very recent experiments have 
generated solitons that oscillate and collide in a harmonic 
trap, in the crossover regime between one dimension (ID) 
and three dimensions (3D) [l7|. As we have shown in a 
recent publication 1^, the decrease of the soliton width 
with the number of atoms (a nonlinear effect) gives rise 
to an improved position resolution, surpassing the shot 
noise threshold, when more and more particles are used 
for the experiment. 

Atom clouds are typically analyzed by taking (pix- 
elized) absorption or phase contrast images, giving access 
to the atomic density profile. If one wants to measure 
some quantity of interest, the resolution is limited by the 
fluctuations in the image due to counting noise on ev- 
ery pixel. We characterize these fluctuations and specify 
how to extract, in an (almost) optimal way, information 
from the noisy data. It turns out that the quantum fluc- 
tuations beyond the Hartree-Fock product state actually 
do not spoil the image resolution, but even improve it 
slightly, although the power law scaling with the atomic 
density (~ n^^' "*) is unchanged. In this paper, we present 
a description of these elements of the theory, providing 



a number of technical ingredients and details that were 
only briefly mentioned in Ref . |18j . 

In Secini we introduce the quantum image of an atomic 
cloud, define the associated Fisher information (FI) and 
recall its connection to optimal parameter estimation, 
starting from the image data. We treat solitons and vor- 
tices in a mean field description with Poissonian counting 
statistics, and we introduce a general Gaussian Ansatz 
for the counting statistics within an image, for which 
the Fisher information can be determined analytically. 
Different scalings with the atomic density are identified 
that range between the shot-noise and the Heisenberg 
limits, depending on the system dimension and the kind 
of nonlinearity. In Sec lIIIi a quantum field theory of 
atomic density correlations is developed within the Bo- 
goliubov approximation. We provide a discussion of the 
role of phonon and zero (or Goldstone) modes, and we 
present detailed calculations for the density fluctuations 
in the image of a dark soliton. The location of the soli- 
ton is, due to the breaking of translational symmetry of 
the problem, itself associated with a Goldstone mode, 
and the corresponding contribution to particle fluctua- 
tions are analyzed. In this section, we also construct a 
nearly optimal protocol for image processing. In Sec lIVI 
we give a brief summary and conclusion. 



II. MEAN FIELD IMAGES 

A. Atom density statistics 

The continuous popularity of cold quantum gas physics 
is due, in part, to the possibility of measuring atomic 
density profiles by fluorescence or absorption imaging. 
Adopting the language of second quantization, the corre- 
sponding observable is the intensity or density operator. 



n(x) = *t(x)#(x). 



(1) 



integrated along a line of sight. (We shall omit this in- 
tegral in the following.) An 'image' thus corresponds 
to a set of measured densities {p(x)|x}. Averaging over 
many images, one gets an estimate for the expectation 
value p(x) — (n(x)). If the detector integrates the den- 
sity signal over some small, but finite area ^px(s) (the 
'pixel' no. s — 1,2,... Af ), one deals with a discrete set 
of operators like 



dxn(x). 



(2) 



.4px(s) 



Correlations between the atomic density in different 
points are related to the field operator in the following 
way, using the bosonic commutation relations 

(n(x)n(x')> = (*^(x)^t(x')^(x')*(x)) 

+ (n(x))J(x-x'). (3) 



This quantity defines the density correlation function 

7'(x, x') = (n(x)n(x')) - (n(x)) (n(x')) (4) 

that will play a key role in this paper. 

Let us illustrate these concepts for the case that the 
quantized field operator ^(x) can be reduced to a single 
mode. This is a common approximation at low temper- 
atures where a macroscopic fraction of atoms condenses 
into a single spatial wave function. The operator ^(x) is 
then replaced by a single annihilation operator oq mul- 
tiplying a classical complex field $(x), and the average 
density is given by 



(n(x)>=A^o|<i>(x)| = 



(5) 



where A^o — i^-o'^o) is the number of atoms in the 'con- 
densate mode' ^(x) (itself normalized to unity). The 
density correlations, essentially the structure factor of 
the system [l^l, are found as 



7'(x,x') ^ {AN^-No)m^rm^'W 
+ 7Vo|<I'(x)|2j(x-x'). 



(6) 



The first term, proportional to the Mandel parame- 
ter [201 , vanishes if the system is in an eigenstate of 
the operator ao, i.e., a coherent state (Poisson statistics). 
The second term describes local fluctuations at the same 
position (for an image: in the same pixel), with a vari- 
ance that is equal to the mean density (the mean atom 
number on the pixel). 

It is intuitively clear that when more than a single 
spatial mode are taken into account in the field operator 
expansion, valuable information about the field's quan- 
tum state is hidden in the correlations of the atomic den- 
sity [2l|. This has been discussed recently for the Mott 
insulator-superfluid transition [22l l23j , and exploited in 
measurements on the strongly correlated Mott insulator 
phase with ultracold ^^Rb atoms released from an opti- 
cal lattice \M\. Excited modes also introduce additional 
fluctuations into an image, however. The competition 
between these two effects will be the central theme of 
this paper. 



B. Parameter estimation 

1. Information measures 

The full joint probability distribution of the atom num- 
bers in every pixel provides the complete 'counting statis- 
tics' of an image. The information content of an image 
is given by the counting statistics via the classical infor- 
mation entropy /, which for a pixelized image is given by 
(e.g., see (H) 



I = - 



Yl P({pJ)l0gp({ps}). 



(7) 



Pi,---Pm 



In ([7]), p{{ps}) is the joint probability measure for the 
occurrence of detection events with p^ atoms detected in 
pixel px(s) (M is the total number of pixels). 

A particular application of this information concept 
appears when we want to estimate a parameter q "hid- 
den" in the image, i.e., the counting statistics is a func- 
tion p{{ps}; q) of the parameter g, referred to as the like- 
lihood function (LF) in statistical estimation theory. The 
hidden parameter can be the fraction of particles in the 
condensate mode, the width of the cloud, or the position 
of a "solitonic" excitation, i.e., the location of a mini- 
mum in the density profile. This will be our example 
throughout the analysis. Information processing theory 
provides an explicit formula for the optimal signal-to- 
noise ratio in measuring q. It can be translated into 
a lower limit on the variance, the Cramer-Rao bound 
(CRB), Var(g) > 1/ F(q ). Here, the Fisher information 
F{q) is given by [H, [2^ 



F{q) 



Yl pi{P'^}'^i) 



d'^\ogp{{ps};q) 



(8) 



In a quantum mechanical framework, the FI has a ge- 
ometrical interpretation as distance (metric), depending 
on a parameter q, in the space of density operators [27|. 
The multiple sum in Eq. ([8]) is difficult to evaluate in gen- 
eral, and we shall focus in this paper on two schemes 
where the calculations are feasible: (i) the single-mode 
approximation for the field operator, assuming Poisso- 
nian atom number counting statistics; (m) a Gaussian 
approximation for the probability measure, where the re- 
sults can be expressed in terms of the average density and 
the density correlations. The latter scheme will be ap- 
plied to a multi-mode field theory within the Bogoliubov 
approximation. 

To estimate the parameter q we require a definite pre- 
scription of how to extract it from the data that fluctuate 
from shot to shot. We shall in particular identify an op- 
timal prescription that permits to saturate the Cramer- 
Rao bound. 



2. Poissonian counting statistics 

The single-mode approximation is well known as the 
mean- field theory for BEC. We assume here that the 
variance of condensate particles is normal, AA'^ = A^o 
(Poisson statistics) so that the density correlations are 
given by the last term in Eq.Q. The mean field theory 
is equivalent to a Hartree product state Ansatz for the 
many-body wave function. As a result, the probability 
measure p{{ps\) factorizes into Poissonian statistics for 
each pixel. The summation in Eq.® can then be per- 
formed analytically. When we take the limit of infinitely 
small pixels, one finds an integral over the spatial coor- 
dinate of the image [28| 



F{q) = 4 I Ax 



d\^{x-q)\ 



(9) 



where we have made the dependence of the average den- 
sity profile [^(a;;^)!^ on the parameter q explicit. (From 
here on, the complex field $ is not normalized to unit 
norm, but its square gives the mean density of the con- 
densate particles.) 

The meaning of this formula can be illustrated by a dis- 
cussion of the optimal signal processing strategy. Given 
the image data p{x;q), we construct, as in Refs. 28', 2!||, 
a linear filter g{x) to provide an estimate for q: 



S{q) 



dxg{x)p{x;q) 



q dxg{x)dqp{x;0), 



(10) 



where g{x) is a local gain function on the pixel at po- 
sition X that can take positive and negative values. In 
the second step we have assumed, without loss of gen- 
erality, that the signal vanishes when the parameter q is 
zero, and we have performed a Taylor expansion of the 
density profile, p{x;q) ~ p{x;0) + qdqp{x;0). The ex- 
pectation value of the signal, S{q), is simply obtained in 
terms of p{x; q) = {n{x)) and its derivative. 

The variance of the signal is found by squaring Eq. (|10p 
(first line) and expressing the average in terms of the 
density correlation function ^ 



AS' 



dxdyg{x)g{y)'P{x,y;q). 



(11) 



Recalling that different pixels are uncorrelated [Eq.® 
reduces to its last term], this variance reduces to 



(12) 



AS' = dxg'{x)pix;0). 



We can now choose the gain function g{x) such that the 
signal-to-noise ratio SNR = 5^(q)/AS'^ is maximized. 
This optimization problem has the following solution, as 
pointed out in Ref . |29[ for a coherent state of light pop- 
ulating a single spatial mode: 

'dmx;q)\' 



gopt{x) 



mx; 



dq 



(13) 



q=0 



Here a is a normalization constant. This gain function 
can also be interpreted as an optimized spatial mode 
(the "noise mode" in the language of Ref. [2^), onto 
which Eq. (fTO|) projects the image. The minimum un- 
certainty Aq for the parameter estimation corresponds to 
an SNR of unity, and this reaches the Cramer- Rao bound 
{Aqf = Var(g) = l/F{q) with F{q) given by Eq.©. If 
the wave function $ were proportional to -s/Aq, we would 

1 /2 

find a shot- noise limited resolution, Aq ex 1/A^o • l-*^^ 
to atom-atom interactions, this limit can be overcome, as 
we shall see. 



3. Gaussian images 

Given the mean atom number per pixel Ps{q) and 
the (pixelized) density correlation matrix Psj (q) with 



[cf. Eq.Q] 



Psiq) 



Psjiq) 



dxp{x;q), 



(14) 



^px(s) 



dx 



dyV{x,y;q), (15) 



^px(s) ^px(i) 



we may make the assumption of a joint Gaussian proba- 
bility distribution. Since a Gaussian is fully characterized 
by its first and second moments the probability measure 
is simply given by |25[ 



p{p\q) 



(2 7r) 



-Af/2 



Vd^t(P) 



cxp 



1 



(p-p)-P-i (p-p) 



(16) 



where the vector p of length M collects the detected atom 
number variables, p collects the corresponding mean val- 
ues, and P is the covariance matrix with elements Psj{q). 
The likelihood function p(p; q) depends on q through 
Ps{q) a,nd Psj{q). 

With the Ansatz (fT6|) for the LF and by replacing the 
sum over discrete particle counts by continuous integrals, 
the FI is given by the following analytical expression |1S] 



92det(P) 




5gdet(P) 
det(P) 



(17) 



E 



92 (P- 



P . + 2fp-M .^M 



If no correlations exist between neighboring pixels, this 
expression reduces to Eq.® for the Poisson case, pro- 
vided one takes both the limit of small pixel size and large 
average atom number per pixel. This is as expected since 
for a large average, Poisson and Gaussian statistics be- 
come similar. We discuss an example taken from Ref. 18 1 
below. 

Correlations between pixels may increase or decrease 
the FI and make the parameter estimate more or less 
precise. A detailed calculation is discussed in Sec lIII Cl 



C. Examples: kinks and vortices 

We evaluate now the Fisher information for measure- 
ments of the position of a solitonic excitation in a Bose- 
condensed gas. We consider first a ID setting where the 
quantum field theory is given by the following nonlin- 
ear evolution equation for the field operator [the Gross- 
Pitaevskii equation (GPE)] [H, [^ 






2 m dx^ 



+ 14, 



3^1/% 



*, 



(18) 



where m is the atom mass, T4xt is an external poten- 
tial, and g the binary interaction strength proportional to 



the s-wave scattering length Og. Within mean- field the- 
ory, the field operator is replaced by the complex order 
parameter $(a;) (the condensate wave function or order 
parameter) that also satisfies Eq. p^ . The mean atom 
density is then p{x) — |$(a;)p, thus normalizing the order 
parameter to the particle density. 

A family of "dark" solitonic excitations of the Bose 
condensate exists for repulsive interactions {g > 0). In 
the homogeneous case (Kxt = 0), it is given by [l3| 

^ix;q) = — liv + \/ c^ - u^ tanh [K(u)(a; - g)] | ,(19) 



which shows a local 
position q = q{t) - 
dip is specified by k{v 



minimum in the density at the 

(?(0) -I- vt. The width of this 

v^ I (? where 



^ — fij \J2mgn is the so-called healing length and n is 
the one-dimensional "background" density. The soliton 
moves with constant velocity v that does not exceed the 
sound velocity c = ^Jng/m. For a soliton at rest, the 
density is strictly zero a,t x = q, and the phase of the con- 
densate wave function increases by tt when crossing this 
point (a "kink soliton"). The kink becomes wider and 
the density dip disappears as the velocity v — > ±c. The 
wavefunction in (jl9[) is normalised such that the atom 
density asymptotically approaches the constant value of 
n (actually, after few healing lengths). Such solitons have 
been created in BEG experiments by shining a laser field 
on one half of the atomic ultracold cloud, which induced 
a relative phase by the AC Stark effect [3l|. 

Topological excitations in two dimensions are vortices 
where the phase of the wave function increases by a mul- 
tiple of 27T when circling around a zero in the density. We 
discuss an example below. 



1. Dark soliton 

Due to the shot noise fluctuations in the detected 
atoms the density dip is washed out and one has to find 
a good estimate for the kink position q. The formula ^ 
for the FI within Poisson statistics can be evaluated ex- 
actly for the dark soliton given by Eq.((T9]). The result 
is [recall that Eq. implies the limit of infinitely small 
pixels] 



47i\/2 V , 
F = ■; arctan 



v/c 



v/c — c/(2 v) 
V2^Kiv) 



(20) 



arctan 



V2e'«( 



8n /„ V , , , 
^2 + -^ ).(.). 



As expected from translation invariance, after integration 
over X, this does not depend on the soliton position q . In 
FigUwe plot the minimum uncertainty Aq — F{q)^^''^ 
of the soliton position in units of the condensate heal- 
ing length, as a function of the soliton velocity v. This 
behaviour makes good sense: when the soliton is almost 
at rest, the dip is very sharp and a precise knowledge of 
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Figure 1: (Colour online) Soliton position uncertainty Aq, in 
units of the healing length ^, versus soliton velocity v (c is 
the speed of sound). The scaled Fisher information is given 
hyF' = Fe- 



the soliton position can be gathered; when the velocity 
approaches c, the dip in the density becomes very shal- 
low, and less information is available from density mea- 
surements. For comparison, the dashed curve shows the 
width 1/k{v) of the density dip: we see that the mea- 
surement precision can largely exceed this value. The 
maximum information is 



v = 



F 



16 n 



(21) 



which implies that, for a soliton at rest, the uncertainty 
in q scales like n~^''^ with the background density n. As 
pointed out in Ref.[i8|, this is a better scaling than the 



usual shot noise limit 



-1/2 



) . We emphasize that this 



enhancement does not require any squeezed or otherwise 
entangled multi-atom state. It simply follows from the 
shorter wavelengths that occur in a BEC matter wave 
interferometer due to the atom-atom interactions. 

The scaling with the density can be understood with 
the following statistical argument: the atom number on 
a pixel of area Ax at position x is given by N{x) — 
p{x) Ax ± yj p{x) Ax, assuming fluctuations at the shot- 
noise level. From error propagation theory the uncer- 
tainty Aq on the soliton position q is given by 



Aq. 



AN{x) 



\(lp{x)/dx\Ax ^^"Ax' 



(22) 



where the healing length ^ sets the scale for the density 
variation around the dip. By processing different data 
points across the relevant region, this uncertainty can 
be reduced by a factor 1/^fM where M « £,/ Ax is the 
number of pixels across the dip, which leads to a precision 
of the order of 



Aq 



oc n 



-3/4 



(23) 



Let us now turn to the Gaussian approximation to the 
likelihood function of an image and discuss the Fisher 
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Figure 2: (Colour online) (a) Rescaled Fisher information as a 
function of soliton position, q, for a Poissonian LF. The solid 
(black) line is the result for pixel size Aa;/^ = 3, the dashed 
(red) line for Aa;/^ = 2, and the dashdot (blue) for Ax/^ = 1. 
(b) Fisher information as a function of soliton position, q, for 
a Poissonian LF, shown for two different soliton velocities and 
fixed Aa;. The solid (black) line is for ti = and the dashed 
(red) line is for v = 0.3 c. 



information (flTj) for the dark soliton. We recover the 
variance and mean of the Poisson distribution by choos- 
ing a covariance matrix with P^j — Var(/9s) ^sj = Ps Ssj, 
with the variance of the number of atoms on the s-th 
pixel equal to its mean value Ps{q)- 

The integral in Eq. p4p with the order parameter given 
in (|19p can be worked out analytically and the derivative 
with respect to q performed. The FI ([T7]), scaled to F' — 
F^^, then becomes 



2 /„2\2 



F' = n'e{l-y /c') 



E 



Je 



9s 



p.(0) 2 ^ p2(0) 



(24) 



where we have put q = without loss of generality and 
gs = {sech^ [k{v){xs + Ax)] — sech^ ['«(u) Xs]} , (25) 



p^(0) = 71 Ax - 2n^'^K{v){tanh[K{v){xs + Ax)] 
— tanh [k{v) Xs]} , 



(26) 



with Xs — sAx, and s an integer. This expression is 
giving those pixels a stronger weight where the density 
profile significantly changes, as could have been expected. 
As a consequence, pixels at the outer limits where the 
atomic density is flat (already a few healing lengths away 



from the soliton position q) do not contribute to the FI. 
We can therefore safely assume that the sohton is located 
well in the interior of a detection window for imaging, and 
even take the limit of an infinitely large window. (For 
a discussion of the limit Ax -^ 0, see Rcf. 18].) That 
is also confirmed in Figl2l^a): we show the dependence 
of the FI on the soliton position q within a pixel. The 
oscillations occur at the pixel size Ax, but for Ax — > 0, 
their amplitude becomes smaller and F' is almost flat. In 
FigEfb), we fix the value of Ax and change the soliton 
velocity v. Also in that case the FI has a maximum for 
q — (soliton at the border between two pixels). 



2. Trapped soliton 

The previous discussion shows that the information 
content is concentrated near the soliton position. This 
suggests that trapped solitons, in an inhomogeneous 
background, behave in a similar way. 

Consider first a soliton at rest in a square potential of 
length 2£. Assuming a single-mode picture with Poisso- 
nian number statistics, the FI becomes 



F'^n^'Y. 



gsjAx) 



^(^-l^^.l), 



(27) 



where the step function ??(x) appears. As expected, only 
those pixels covering the confining box contribute to the 
information. In the limit ^ — > oo, we recover the result of 
an infinite Bose gas given by the first sum of (|M|) with 
V = Q. The same discussion can be applied to the Gaus- 
sian approximation for the complete counting statistics. 
In a harmonic trap, the order parameter can be ap- 
proximately written as [32| 



^{x;q) ~ $bg(a;) tanh 



VV2ec 



(28) 



where $bg(a;) is the Thomas- Fermi (TF) solution for the 
trap ground state [33[, ^o is the healing length for the 
background density at the soliton position (it depends on 
the density uq — |$bg(9)P [SJ])- Let us focus on solitons 
at rest, close to the centre of the trap, q ~ 0, and much 
smaller than the size Rtf of the background field. The 
FI ^ then turns out to be 



-Ptrap(g) 



16 no 



(29) 



which is very similar to the homogeneous case, keeping 
the density no the same. 

This result can be written in different ways by express- 
ing the central density no in terms of other parameters. 
Taking, for example, the total number of atoms, we find 



Ft 



trap 



14.1 



iVo 



(30) 



where Ux is the (single-particle) ground state size in the 
harmonic trap. This leads to the usual shot-noise scaling 

— 1/2 

Nq for the soliton position. Most of the atoms, how- 
ever, are not needed for the image since the information 
content is concentrated near the soliton minimum. In this 
respect, the trapped situation is favorable compared to 
the homogeneous case since the central density is higher 
and the healing length is smaller in the center than in 
the condensate wings. This can be made quantitative by 
comparing the harmonic trap to a homogeneous sample 
of the same size and same total number of atoms. The 
ratio of the Fisher Information between the two cases is 
then 



Fho 



trap 



0.11, 



(31) 



where corrections of order ^/Rtf are neglected, consis- 
tent with the TF approximation. It should be noted that 
the size of the system in both cases scales in a different 
way with the atom number. 



3. Solitons with three-body interactions 

The GPE relies on the two-body pseudopotential 
g6{x), which describes the interaction between the parti- 
cles. One can also consider low-dimensional Bose super- 
fluids in a strong-coupling regime with three-body colli- 
sions. We analyze here the solitonic solution in a quasi 
ID setup governed by the equation (homogeneous case) 






2 TO dx'^ 



im' 



$ 



(32) 



with the universal coupling constant 7 = (fi.7r)^/(2rn). 
This system has been investigated in detail in Refs.f35|, 
l3a |. There exists an analytical soliton solution for the 
order parameter ^{x;q) = ^/nf{x)e^'^'^^\ with 



f{x-q) 



i[l-{v/cf] 



^/l + i{v/c^ cosh[n' [x — q)] 



27rnVl - {v/cf 



(33) 
(34) 



This kind of solitons can also appear in a degenerate elec- 
tron plasma [37| , since the equation for the order param- 
eter describing the electron density is almost the same. 
As opposed to Eq. (fT9)l . the soliton width 1/k is scaling 
here r^ \jn with the density n, and we therefore expect 
an improvement in the soliton kink estimation. Indeed, 
using the Poisson formula ^ for the Fisher information 
we obtain for a soliton at rest 



F = 2%/37rln(2-|-\/3)i 



(35) 



that can be compared to Eq. ([2T|l . With such topological 
excitations we thus reach a 1/n scaling as in the Heisen- 
berg limit, without any entanglement or squeezing. 



4- Vortex line 

Finally, we consider a two-dimensional situation with 
a vortex line. The vortex solution of the 3D GPE is de- 
scribed by a core region where the density goes to zero. 
Vortices are not stable solutions, however, and only in a 
frame rotating with a high angular velocity they corre- 
spond to local or global energy minima [l3]. Vortices in 
a rotating BEG with respect to the lab frame were ob- 
served in several experiments, as for example in Ref.[38|. 

We assume the simple case of a Bose gas with uniform 
confinement along the z axis (length L) and homoge- 
neous in the xy-plane. The solution of the GPE for a gas 
rotating around the z axis is given by [19[ 



$(r;g) = ^f^ f{x,y,z\q) 



X — q + iy 
x-q-iy 



,/2 



(36) 



where n is the two-dimensional density, / is a real func- 
tion, {q, 0) is the position of the vortex line in the xy- 
planc, and s is an integer (the 'winding number' or 'topo- 
logical charge'). The order parameter ([5S|) is an eigen- 
function of the angular momentum with eigenvalue hs. 
Since we are interested in the calculation of the FI ^, 
which depends only on the absolute value |$|, the angu- 
lar dependence is not relevant in the present setting. In 
the following we specialise to the case of a singly charged 
vortex s = ±1. The function / can then be well approx- 
imated by [30l | 



f{x,y,z;q) 



^2e + {x-qy+y^ 



(37) 



We generalize Eq.® for the FI to this two-dimensional 
setting and obtain 



F ^ 4 dz dxdy 



d\<i>{r,ip,z; 
dq 



nn, 



(38) 



where the integration is easy to perform in cylindrical co- 
ordinates. It is interesting to note that for a vortex line, 
the FI depends on the (two-dimensional) density n lin- 
early. Hence, the fact that the vortex core scales with the 
healing length does not improve the shot-noise limit for 
its position detection. This can be understood from the 
simple argument of Eq. (|23p because the relevant number 
of atoms scales with n^^ ~ h /(2mg2D) which is only 
weakly dependent on the density: the dimensionality of 
the topological excitation plays a crucial role. 



5. Dimensionality crossover 

The last example has shown that the dimension of 
the topological excitation is very important. Coming 
back to the quasi-lD situation, we discuss here briefly 
the crossover of the nonlinear terms in the GPE as one 



changes the confinement transverse to the long trap axis. 
The ID solitons described by Ea.(fT5|) are difficult to re- 
alize in an experiment because one needs a high ratio be- 
tween radial and axial trap frequency. When transverse 
effects are taken into account with more accuracy, the 
GPE is modified into a ID non-polynomial Schrodingcr 
equation [33, EO, Ell, E^]. This changes, for example, 
the oscillation frequency of a dark soliton in a harmonic 
trap [4l| . In the limits of weak and strong coupling, the 
nonlinearity in the effective GPE becomes polynomial, 
and analytical results for solitonic solutions (bright and 
dark) can be given [39J. By inspecting these dark soli- 
tons, it is easy to show that the Fisher information scales, 
in the strong coupling regime, as n^ with the background 
density n, identical to the scaling found in Sec lIIGSl 
This beats the shot noise scaling {F ex n) and illustrates 
that also with the help of two-body interactions, one can 
reach the scaling of the Heisenberg limit in a strongly 
interacting system. 



III. QUANTUM FIELD THEORY FOR IMAGES 



Motivation 



Mean field theories are ubiquitous in physics and play 
an important role in the explanation of many phenomena 
in condensed matter physics [43, E3, e.g. superconduc- 
tivity in metals, and even in high energy physics as for 
example the quark condensate in the so called instan- 
ton ensemble in QCD (Quantum Chromodynamics) [45| . 
The mean field approach describes successfully key fea- 
tures of quantum degenerate dilute gases [l^ [3y|, for 
example the static, dynamic, and thermodynamic prop- 
erties of trapped Bose-Einstein condensates [i^, and it 
has been well confirmed in a number of experiments (see, 
e.g., Ref . [47| ) . The mean field approximation applied to 
an ultracold Bose gas leads to the Gross-Pitaevski equa- 
tion ([THj). The stationary solutions of the GPE repre- 
sent the macro scopically occupied spatial mode functions 
when the temperature of the trapped gas is well below 
the critical transition temperature. This concept can be 
formally translated, in the C/(l) symmetry broken ap- 
proach, in the following way: the total matter field is 
split into ^(x) — $(a;) + 5'^{x), where 5'^{x) represents 
quantum fiuctuations around the Gross-Pitaevskii solu- 
tion ^{x). The mean field approach essentially neglects 
the fluctuation field. 

Mean field methods, however, assess only the average 
atomic density and do not provide information about 
quantum noise correlations, which are important in the 
understanding of quantum phase transitions in ultracold 
atoms, e.g.. anti- ferromagnetic structures or charge den- 
sity waves [48 1 . 
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B. Bogoliubov approximation 

1. Phonon and zero modes 

Even at zero temperature, the particle density shows 
fluctuations that lead to the phenomenon of "quantum 
depletion" [T^, S, [50|- This can be described by lin- 
earizing the total many-body problem around the mean 
field solution. In other words, one expands the second 
quantised Hamiltonian up to quadratic terms in the fluc- 
tuation field (5^. The dynamics of the latter is then gen- 
erated by the Bogoliubov-de Gennes linear operator 



L = 



n 






5 $2 



(39) 



where we have introduced the Gross-Pitaevskii Hamilto- 
nian 



n^-:——^ + V,.tix)+gmx)f-fi, (40) 



and fi is the chemical potential. 

The mode expansion of the fiuctuation operator 6"^ 
is complicated by the fact that C has eigenvectors with 
eigenvalue zero. These arise from continuous symmetries 
that leave invariant the energy of the Bose-Einstein or- 
der parameter <I> [4^. A well-known example is the U{1) 
global phase invariance of the GPE (|40|) , another one is 
the translation in space of the soliton solution (fT9|) for a 
homogeneous system. By choosing a dark soliton solu- 
tion with a definite phase 6 and a given position q, one 
has spontaneously broken these symmetries. According 
to the non-relativistic Goldstone theorem [4J, IMI, this 
spontaneous symmetry breaking is associated with gap- 
less excitation modes, the Goldstone modes, for any sys- 
tem of particles with finite range interactions. The Gold- 
stone modes produce quantum fiuctuations of the phase 
and position of the soliton order parameter, even at zero 
temperature. The global phase is conjugate to the num- 
ber of particles, and the phase fiuctuations can indeed 
be interpreted as a consequence of the variations of the 
chemical potential with particle number [331] . 

The eigenvalues Ek of the Bogoliubov-de Gennes op- 
erator p9p correspond to eigenvectors (uk,Vk) that we 
denote "phonon modes" . They are normalized and or- 
thogonal according to 

Skp^ dx {ul{x)up{x) - vl{x)vp{x)} (41) 



and are accompanied by partner modes (wfc*,Wfc*) with 
negative norm and eigenvalue —E^- This construction 
fails for the so-called zero modes (uq,, Wq) that are in the 
kernel of C [50, [52, Hj]. For the global U{1) symmetry, 
the zero mode is ($, — $*) that is generated by applying 
-id/d0 to ($e*^ $*e-*^). The norm gU of this mode is 
zero. We need partner or adjoint modes {u'^, v^) to sat- 
urate the completeness relation in the space of fiuctuation 



fields. They can be constructed by solving the equation 
C{u^,v^'^) = ma{ua,Va). The generalized eigenvalue 
nia is called the "effective mass". One gets zero modes 
of C^ and a completeness relation in the form [53| 



^i^ ~y) ^ "^{Mx) uliy) - vlix) Wfc(y)} 
fe 

+ J2{uo.{x)K\y)]*~vf{x)v*Jy)}, 



(42) 



where a enumerates all broken symmetries and the ad- 
joint modes are normalized such that 



6^0 = / da; {ul^*ix)uisix) - v%^*ix)vpix)} 



(43) 



In our one-dimensional soliton system, we expect two 
Goldstone modes associated with the global soliton phase 
9 and its position q. The corresponding quantum fiuctu- 
ations will prove important for our analysis. We give the 
phonon and zero modes for this geometry in Appendix lAl 
The fiuctuation operator S^{x) around the Gross- 
Pitaevskii solution ^{x) is expanded as 



S'i'ix) =J2{bkUk{x) +bivlix)} 

k 
+ ^{Pc.<'^(x) - i gaUc(x)}. 



(44) 



The operators are constructed to implement the commu- 
tation relation [5^! [x) , 5^!^ {y)] = S{x — y): hk, b\ are 
bosonic annihilation and creation operators, [&fc,fol] = 
(5ap, and for all broken symmetries in the problem, 
the "position" and "momentum" operators Qq, Pa are 
canonically conjugate, [QonPfi] = i^ap- The complete- 
ness and orthogonality relations ((42l [43|) are compatible 
with the construction 

Q« = I lAx[uf*{x)5i>{x)~v^^*{x)6i!\x)^iAh) 

The operator Pa can be identified with the generator of 
the symmetry transformation behind the corresponding 
zero mode. 



2. Mean density and correlations 

We now provide a general framework for the mean 
atomic density and the density correlations within Bo- 
goliubov theory. Going beyond Ref. [l^l , we consider the 
case of phonon modes at finite temperature. The field ex- 
pansion (|44|) is useful for the computation of expectation 
values because different modes are not correlated. 



Evaluating the average density, we get 

k 

+ J2{bih)\uk{x,q)f+Z{x), (47) 



that appear Hke "mode functions" for the zero mode 
operators. We have simpHfied Eq. (P^ by noting that 
Im [m* (a;) M„'^(a;)] = for the zero mode functions given 
in Appendix [Xl To proceed, we have to specify the quan- 
tum state of the sohton displacement modes. 



with the occupation number (h\.bk) = \/{e^^'' — 1). For 
thermal states, (&g) — 0. The so-called quantum deple- 
tion is related to the finite contribution ~ |wfcP of phonon 
modes even at zero temperature: this is a direct manifes- 
tation of quantum density fluctuations. Finite tempera- 
ture adds thermal contributions to the phonon modes. 
The zero mode contribution is obtained as 



Z{x) = Y. {l<'(^)l' ^Pl) + \^o.{x)\'' {QD 



(48) 



-Re K{x) u:f{x)] - Im K{x) <'^(x)]({P„, Q4)} , 

where {•,•} denotes the anticommutator. For the global 
C/(l) symmetry, the relevant operator averages are: 
(P,2) = No, {QD = 1/(4 TVo), and ({Pe, Qs\) = 0, they 
follow from the assumption that the condensate mode is 
in a coherent state with A^o particles on average. The 
state for the Goldstone mode associated with soliton dis- 
placement is discussed in Sec lIII B"3l below. For the mo- 
ment, we only assumed that the operators Qa, Pa av- 
erage to zero which is plausible since they appear only 
quadratically in the Hamiltonian. 

When quantum fluctuations around the mean field are 
taken into account, different spatial locations can become 
correlated beyond the level of Eq. ([6]) because they probe 
the same delocalized excitation modes. We find the den- 
sity correlation function V{x, y; q) [Eq.(|3])] by a straight- 
forward expansion of the four-point field correlations to 
second order in the fluctuation operator 5'^ (this is con- 
sistent with the Bogoliubov approximation) . Thus we get 
for phonons in a thermal state 

V{x, y; q) ~ ([$*(a;)(5*(a;) + 5'^\x)^{x)] 

x[$*(y)^*(y) + ,5*t(y)$(y)]) 

= E{(i + (^I^^-))^(^)/fe(2^) 

k 

+ {blik)h{y)f'k{x)] 

(49) 



where we introduced the functions 

fk{x) = •i>*{x)uk{x)+Vk{x)'^{x) 
fl{x) = ^*{x)vl{x)+ul{x)^{x) 

for the phonons, and the abbreviations 

7?„(x) = '^*{X)U^^{X)+U^^*{xmx) 
iPa{x) = i^*{x)Ua{x) - iu*^{x)^{x) 



(50) 



(51) 



3. Quantum statistics of soliton position 

In the second-quantized many-body theory, the soli- 
ton position is described by an operator Qq whose fluc- 
tuations "fiH" the density dip, as discussed by Dziar- 
maga [S^, l5J] . In a homogeneous geometry, the dynam- 
ics of Qq is similar to a free particle, the correspond- 
ing Hamiltonian in the Bogoliubov approximation being 
Pq/2mq. The stationary states of this Hamiltonian are 
momentum eigenstates where the soliton position is max- 
imally uncertain. However, the effective mass occurring 
here is negative, iriq = —'im{n/n). This mimicks the 
negative kinetic energy of the classical (non-quantized) 
moving soliton |55l . |56[ . A similar phenomenon occurs in 
a harmonic trap where the soliton displacement mode has 
a negative frequency (the mode is called anomalous) [33| . 
These degrees of freedom are therefore thcrmodynam- 
ically unstable, and we cannot use thermal statistics. 
Dziarmaga in Ref.[32] has suggested an alternative ap- 
proach to specify the quantum state of this degree of 
freedom, to be used as initial condition for the subse- 
quent dynamics. The idea is to "pin" the soliton to the 
position q by minimizing the density ^^((7)^(5), which is 
a quadratic form in Pq and Qq [see Eg. pS)) aX x — q]. As 
shown in Ref. [l8| , one finds in this way a Gaussian state 
similar to the ground state of a harmonic oscillator with 
respect to the "Hamiltonian" (in quotes since it has not 
the dimensions of energy) : 



h{Pq,Qq) 






riK^Ql 



K 

4' 



(52) 



where the mode functions of Appendix [X\ Eqs. llATjIASp 
have been used. The ground state of this Hamiltonian 
gives (h) = and can also be written as the vacuum 
state corresponding to the annihilation operator 



(53) 



With the help of bq and its conjugate operator, we can 
write the expansion of d'^{x) in the same form as for the 
phonon modes. 

We consider here two classes of states that general- 
ize this ground state: "squeezed" states and "thermal" 
states. The squeezed state depends on the positive pa- 
rameter C: it is defined as a Gaussian state with quadra- 
ture variances 



{P^)C 



2nK 



{Qlk 






(54) 
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We note that these states give average density profiles 
that resemble a partially filled dark soliton, as illustrated 
in FiglHl The deviations from the "optimal case" are con- 
trolled by the parameters C or r. The quantum states 
constructed here are not stationary states of the Bogoli- 
ubov Hamiltonian, however, and will evolve in time, as 
discussed in Ref . S2I ■ 



Discussion of results 



Figure 3; (Colour online) Density profiles in the notch: the 
solid (black) line corresponds to the squeezing parameter (^ = 
1 whereas the dashdot (red) to ^ = 100, while the dashed 
(blue) line corresponds to the thermal state with parameter 
T = 5/^. Other parameters are: ^ = 10^, n^ = 100, and 140 
phonon modes have been taken into account. 



where C = 1 corresponds to the ground state. 
For this state, we get an average density {h)(; = 
(k/8) (C"^ + C - 2) > 0. The density at x ^ q van- 
ishes for C = 1, and FigE] actually shows the result of 
Ea. (|47p . including contributions from quantum noise, for 
that case. 

The thermal state is defined by analogy to the canon- 
ical ensemble as a state that maximizes entropy at a 
given mean value of h. This mean value is given by 
{h)r — |K/(exp[K/(2r)] — 1) and it is controlled by a 
parameter r with the dimension density. The quadra- 
tures are 



(p: 



2nK coth[K/(4r)], 



(Ql) 



3th[K/(4r)] 



8nK 



In Fig[3]we show density profiles corresponding to dif- 
ferent values of the squeezing parameter C = li 100 and 
for the thermal parameter r — 5/^. The plot confirms 
that the choice C = 1 (or when r ^ 0+ for a "ther- 
mal" q-statc) gives the state with minimum depletion 
that best resembles a condensate with a perfect soliton. 
For r = 5/^ the minimum in the notch is not zero (as for 
C 7^ 1), but the density at the edges becomes higher than 
the density of the squeezed quantum state, refiecting the 
"thermal" and the long wave nature of the zero q-mode 
state. 

For what concerns the density correlations, it is con- 
venient to bring them in a form where the (5-correlated 
term in Eq.Q that appears due to normal ordering, is 
subtracted. To this effect, we use the completeness re- 
lation (|1^ to rewrite the sum over Uk{x)u*f.{y) in the 
correlation function (l49l) as 



V{x, v) = $(a;)$(2/) [5{x - y) + J{x, y)] , (56) 



(55) where the function J^{x,y) is found as 



Jix,y) =Re< ^[2vk{x)vl{y) + Ukix)vl{y)+Vk{x)ul{y)] \ 

+ 4{{P^)uf{x)uf*{y) + {Ql)u,{x)u;{y)} - J2{u^{x) K^iy)]* + ul^{x) <(y)}. 



(57) 



r 



In Fig|4]we show the behavior of J^{x, y) for two different 
values of the squeezing parameter. Between C = 1 (left) 
and C, — 100 (right), the overall magnitude of J'{x,y) 
changes significantly, demonstrating a strong influence of 
the soliton displacement mode. For a large uncertainty in 
Qq {( = 100), the function ^{x, y) is positive everywhere 
and concentrated near x k, y k, q. The choice C = 1 leads 
to negative values (anticorrelations) along the diagonal 
X = y away from the point x — y ^ q. 

The calculation of the Fisher information is done nu- 



merically by using the phonon and zero modes of Ap- 
pendix [XI normalizing them in a box of length 2i. In 
FigEJa), we show the Fisher information for different ap- 
proximations: a single-mode BEC with Poisson statistics 
[dotted line: homogeneous system, solid line: in a finite 
box, Eq. ([27|) ] is compared to a multi-mode calculation 
with a Gaussian counting statistics (dashed, dot-dashed, 
and thin solid lines) . The correlation matrix gives contri- 
butions to the FI via the derivatives of P""*" and det(P) 
with respect to the soliton position [see Eq.pT))]. We 
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Figure 4: Plot of the function J{x,y) defined in Eq. (|57|) . Parameters: (. = 10^, n^ — 100, and 140 phonon modes have been 
taken into account. On tire left, squeezing parameter C = 1 (optimal) and on the right C, = 100. 



have chosen the 'optimal dark soliton' C = 1.0 (dashed 
line) and a slightly filled (squeezed) one, C, = 1.5 (dash- 
dotted line). The thin solid line corresponds to the FI 
for a thermal state with a parameter r = 0.2/^. 

First, it is interesting to note that one gets more infor- 
mation in a box than in the homogeneous case. This is 
a finite-size effect which enhances the impact of "missing 
atoms" in the soliton center. Second, the exact result 
for the Gaussian multi-mode theory ([TT]) shows that cor- 
relations, including the zero modes, increase the level of 
information that we can extract, beyond mean-field the- 
ory (single- mode approximation). Moreover, the dashed 
curve for C — 1 shows that minimizing the quantum "fill- 
ing" of the notch provides the highest information. The 
divergence of F' at small pixel size is due to the failure of 
the Gaussian approximation which becomes unphysical, 
as the average atom counts per pixel drop below unity, 
seeRef.'ls]. 

In FiglSJb) we show the scaled Fisher information ver- 
sus the scaled density n. Both mean-field and multi-mode 
theory give a linear scaling with n which translates into 
the same power law n~^''^ for the sensitivity of the soliton 
position. 

In order to understand why the inclusion of noise with 
the Bogoliubov description provides an even better reso- 
lution than the mean field, we use the formalism of signal 
processing introduced in Sec lIIB 2l The signal-to-noise 
ratio gives an assessment of the amount of information 
that can be extracted from a given statistical estimation 
strategy. By using the completeness relation ([1^ and 
the result ([57]) it can be easily shown that the noise AS"^ 
[Eq.inD] splits into 



A52 = A^^^F 



^Sly, - A5^ 



(58) 



for any gain function g{x). The three terms here are 
the contributions of the mean field, the phonon contri- 
bution, and the sum over the Goldstone (or zero) modes. 
The Goldstone contribution is negative and reduces the 
noise even below the phonon level. This is related to 
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Figure 5; (Colour online) (a) Rescaled Fisher information 
F' — F^ versus pixel size for different regimes: solid (black) 
line for a Poissonian counting statistics in a box and the dot- 
ted (green) line in the uniform Bose gas limit. The other lines 
take into account quantum and thermal density fluctuations, 
using a Gaussian approximation to the counting statistics. 
They differ in the quantum state of the Goldstone mode as- 
sociated with the soliton displacement. Dashed (red): state 
with zero density at soliton position ("squeezing" parame- 
ter ^ = 1, effective "temperature" r = 0); dashdot (blue): 
squeezed state C, = 1.5; thin solid (magenta): thermal state 
r = 0.2 /^. (b) Rescaled Fisher information versus linear den- 
sity: dashed (black) line corresponds to Eq. (|21|l . and the dash- 
dot (red) line is the information extracted from the signal-to- 
noise ratio for a gain function g^° (see text), and the solid 
(blue) line is the Fisher information within Gaussian and Bo- 
goliubov approximations for a finite pixel size Aa; = 0.7 5, and 
(^ = 1. In both pictures the numerical simulations were made 
with a box length 2£ = 20^, and by considering 140 phonon 
modes. 



the behaviour of the function J(x^y) in ([57| for C = 1; 
FigHfleft): its negative values are larger in magnitude 
that its positive ones. Recall that for the density correla- 
tions ([SG]) . J(x, y) is multiphed by the product <l>(x)<I'(y) 
that is zero for x = y = q and increases to a positive con- 
stant n for X — y ^ q after a few healing lengths. This 
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further enhances the negative contributions of J'{x,y), 
and reduces the noise AS^ relative to the mean-field ap- 
proach. 

The above discussion clarifies how the Cramer-Rao 
bound can be smaller for multi-mode fields compared to 
the mean field description. But can one also identify a 
strategy to reach the CRB? In the Poissonian case, with- 
out correlations, the SNR with the optimal gain func- 
tion ^opt (x) [Eq. P^ ] reaches the sensitivity given by the 
Fisher information ([5]), and this happens when the vari- 
ance A^*^ ex S{q). We shall proceed in the same way with 
the multi-mode theory, but work with finite pixel areas 
to avoid the breakdown of the Gaussian approximation. 
In that case signal and noise are given by 



S{q)^g-^Ax, 

AS'2 = g•PgAx^ 



(59) 



where the vector g (length M) represents the values of 
the gain function g{xs) on the pixels. These expressions 
arc the discrete versions of Eqs. pUlfTTj) . 

We can achieve AS^ ex S{q) by solving the linear sys- 
tem P g = dqp, giving an optimal gain function g^°. It is 
then easy to check that the last term in Eq. (fT7)l becomes 



^^ '"' dq dq dq 



P-'^^S{q). (60) 



factor n"^/'' the scaling of the classical shot noise limit 
1/y/n. We emphasize that this limit is reached without 
the need of any entanglement or squeezing of the system 
state. Even more favorable scalings are found for strongly 
interacting systems where the nonlinearity appears with 
a different exponent in the GPE. In optical and atomic 
interferometry, shorter wavelengths provide a better res- 
olution of phase changes, and we can explain our result 
as a consequence of the high wave number content, i.e., 
the steepness, of the soliton dip. A similar improvement 
could also be obtained by simply applying fast counter- 
propagating beams of atoms. We believe, however, that 
the stability properties of solitons and the fact that there 
is no or only little net particle current makes this sys- 
tem suitable for interferometric investigations over longer 
time scales compared to a thermal beam. 

We have investigated the influence of quantum fluctu- 
ations with the help of Bogoliubov theory in the weakly 
interacting case. The scaling law (n^^/^) remains stable, 
but the prefactor is different. The Goldstone modes that 
are associated to spontaneously broken symmetries, in 
fact increase the sensitivity of the measurements. While 
signal processing theory provides a theoretical limit for 
the sensitivity (the Cramer-Rao bound) , it generally does 
not provide a method to achieve this limit. In our ap- 
proach, we were able to find a nearly optimal filtering 
function, and showed with a signal-to-noise ratio analy- 
sis that it almost reaches the Cramer-Rao bound. 



We were not able to find a similar proportionality for 
the term Ti[{d'^P^^ /dq'^)P]. A numerical analysis, illus- 
trated in FigjSJb), shows that we are nevertheless very 
close to the CRB with the gain function g^°. 

We emphasize that this argument only requires the 
first and second moments of the counting statistics to 
optimize the SNR. It does not directly rely on the Fisher 
information which is only available within the approxi- 
mation that the entire counting statistics is Gaussian. 
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IV. CONCLUSION 

We have analyzed in this paper the quantum fluctu- 
ations in the counting statistics of an atomic density 
image. We used Bogoliubov theory, and we observed 
that zero modes, due to the breaking of the U{1) and 
translation symmetries, makes quantitatively a signifi- 
cant change in the information content of the counting 
distributions. We applied our theory to identify the ulti- 
mate information theoretical limits for position measure- 
ments of dark solitons and other topological excitations 
in a quasi one-dimensional and two-dimensional Bose- 
Einstcin condensate. 

In the case of a pure condensate, where all particles oc- 
cupy the same quantum state described by the GPE, the 
best estimation of the soliton position has an uncertainty 
that scales, for a weakly interacting system, as n~^''^ 
with the linear background density n. This beats by a 



Appendix A: EXCITATION MODES FOR A 
DARK SOLITON 



Here we sketch the derivation of the Bogoliubov eigen- 
modes for a system with a soliton in a box of length 
2£. The analysis complements results given in Ref . |54l| . 
We are looking for the eigenstates of the Bogoliubov- 
de Gennes operator ([55)1 with V;2xt(a;) — 0. The back- 
ground or condensate field is given by the wave function 
$ = ^(x; q) of Eq.((T9]). We focus here on the stationary 
case (soliton velocity f = 0). Simple and compact ex- 
pressions for the Bogoliubov eigenmodes can be found in 
the limit £ ^ £,, when the soliton is located well within 
the quantization box. The integral of |<I>(a;;g)p over the 
box gives the number of condensed atoms Nq: 



No = 2£n - 2n/K 



(Al) 
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where the negative correction describes the atoms "miss- 
ing" in the soliton notch. We have set k = k{v = 0) = 
l/\/2f. For the sake of simpUcity in the notation, we 
shaU hereafter drop the parametric dependence on q in 
the order parameter $, the phonon modes Uk, Vk, and 
the zero modes. 



1. Phonon modes 

Following the approach in Ref. [SJ, [53| , we find that the 
modes of £ in a box with periodic boundary conditions 
can be written as 



Uk{x) 
Vk{x) 



Ml, e' 



+ f3t 



k 

— sech [k (x — q)] 

K 



h i tanhfK (x — q)] 

2k 



(A2) 



where the upper [lower] sign applies to Uk [wfc], respec- 
tively. Here, M^ is a normalisation constant [given in 
Eq.dMI], 




Figure 6: (Colour online) First three Bogoliubov mode func- 
tions ffe(x) for g = 0, normalized in a box of length 2^ = 20^ 
(periodic boundary conditions). Thick (thin) lines represent 
the real (imaginary) part oivk{x), respectively. Solid (black) 
line: fci = 0.5379 (n/E), dashed (red) line: fe = 1.6093 (n/e), 
dashdot (blue) line fca = 2.6704 (tt/^). These fc- values can be 
compared to a free particle where kj = 0, 1, 2 (n/l). 



Pi 



± 



2 Eh 



(A3) 



c is the speed of sound, and the phonon energy is given 
by 



Ek = hc\k\\l + 



fc2 

4^ 



(A4) 



This energy is the same as on a homogeneous background 
condensate. 

The presence of the soliton becomes manifest in the 
total phase shift of a phonon passing from the left to the 
right end of the box. It is given by 2k i + A(p(/c), where 
A(^(fc) = 2 arctan(2 n/k) is due to the interaction with 
the density notch. The quantisation condition for the 
wave number k is thus: 



2kje + A(p{kj) = 2ttj 



J=±l,±2, 



(A5) 



[The case j ^ is excluded because \2k£ + A(p(fc)| > n.] 
In Sec lIII B1 sums over k are understood as running over 
this discrete set of wavenumbcrs. 

From the normalization condition (|4ip for the phonon 
modes, we obtain the normalization constant Mk 



Mk = — 



2k V 2e 



Ik 



2k2 



-1 



-1/2 



(A6) 



In Fig[6] we show the first three cigcnfunctions Vj with 
J = 1, 2, 3, for q = 0. 



2. Zero modes 

For our quantum degenerate Bose gas the Goldstone 
modes originate from the breaking of the global U{1) 



phase symmetry and translational symmetry by assign- 
ing a phase 9 and a value of the displacement q to the 
order parameter $ respectively. In order to determine 
the associated mode functions, we simply differentiate 
the order parameter $ with respect to its global phase 9 
or the parameter q, as described in Sec lIII Bl and Ref. |54l| . 
We obtain 



ue{x) = <^{x) 

Uq{x) — — iK\/nsech [K,{x — q)\ 



(A7) 



with Voi{x) = — u* (a;). We have chosen here 6* = 0, i.e., a 
real- valued $(a:), otherwise the phase of Uq{x) would be 
different. With respect to the Bogoliubov scalar product 
Eq. (|4Tl) . the modes (|A7|) have zero norm and are mu- 
tually orthogonal. The zero modes are accompanied by 
adjoint vectors that satisfy Cu'^ (x u^ and are found as 



uT{x) 



uf{x) 



2{Nqk + n) 



[^{x) + ixUq{x)] , 



4xA^ 



(A8) 



They have zero norm because v^{x) 



,ad* 



(x). The 



adjoint modes are mutually orthogonal with respect to 
Ea. (|4T)) which explains the presence of the second term 
in Eq. llASp . For the normalization and this extra term, 
we consider the limit k£ ^ 1 and neglect exponentially 
small corrections. 
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